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ABSTRACT 


A global, six-level primitive equation model was numeri- 
cally integrated to examine the interaction of a planetary 
wave with a growing mid-latitude baroclinic wave. The 
purpose of this study was to determine the contribution of 
the baroclinic wave to the initiation of a monsoon surge. 
Three experiments were run. Experiment I depicted the 
effect of the planetary wave on the southward migration of 
cold air. The cyclogenesis of a mid-latitude baroclinic 
wave was shown in Experiment II. Experiment III combined 
the growing baroclinic wave with the stable planetary wave 
in an attempt to simulate a cold surge. Characteristics of 
a monsoon surge were present in the numerical simulation but 
not to the extent or detail of observed surges. The results 
indicated the important role of planetary wave advections 


in the formation of cold Surges. 
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I. INTRODUCTION 


The Asian winter monsoon has been an area of very active 
research over the past few years. Winter monsoons are 
dramatic manifestations of land-sea-air interaction. These 
phenomena consist of a mixture of all scales of motions, and 
the various interactions are still not well understood. 
Although the fundamental structure of the monsoon has 
frequently been observed, the dynamical mechanisms are 
still unknown. Recent areas of investigation include the 
role of the winter monsoon in mid-latitude tropical inter- 
action, monsoonal variation and the forcing mechanisms for 
small scale monsoonal variations. 

The monsoon is a three-dimensional, planetary scale wind 
regime which exhibits a strong seasonal dependence. It 
exhibits two distinct seasonal low-level wind flow patterns-- 
a winter outflow from a cold continental anticyclone and a 
summer inflow into a continental heat low. The winter anti- 
cyclone and the summer heat low result from land-sea tempera- 
ture differences and the varying position of the sun. The 
Shapes of the continents and their topographies also play 
a role in producing the considerable regional and temporal 
variability of the monsoons. 

Though centered far to the north, the Siberian anti- 


cyclone dominates the Southeast Asian winter (Ramage, 1971). 
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The low heat capacity of the land relative to the open ocean 
to the east and south ensures the development of the cold 
anticyclone with surface flow from land to ocean. A com- 
pensating flow from ocean to land occurs in the upper 
troposphere. The Himalayas block the southward movement of 
the extremely cold air from the Siberian anticyclone. As the 
anticyclone becomes established Over Siberia, a SUrLace 
cyclone develops over Indonesia and Malaysia as a major 
planetary scale convective area is shifted eastward from its 
Summer poSition near India (Baker, 1983). 

Chang and Lau (1980) analyzed the monthly mean velocity 
potential field at 200 mb for December 1974 and showed a 
broad region of divergence over Indonesia. This broad area 
of divergence iS indicative of deep rising motions and large 
scale overturning in the atmosphere. A region of maximuzin 
convergence over northern China, in the vicinity of the 
Siberian anticyclone, was also noted. 

A strong Hadley cell dominates the winter tropospheric 
circulation of this region. This thermally direct circu- 
lation involves the descent of cold air over the Siberian 
anticyclone, the southward flow of air at the surface, the 
ascent of warm air over Indonesia and a northward return 
flow at 200 mb. 

Although the large scale features of the monsoon are 
identifiable from year to year, there is large interannual 


variability in strength and location. Velocity, temperature 


itil 


and pressure fields vary on short time scales. One such 
variation is the monsoon surge, also known as the "cold 
surge." These events occur at intervals of several days 

to two weeks during the winter monsoon (Lim and Chang, 1981). 

Observational studies by Danielson and Ho (1969), Chang 
et al. (1979), Murakami (1979), Chang and@iau (U930) 1927 
and Chang et al. (1983) discuss Elewe namatere tire clomid 
effects of the monsoon surge at the synoptic and planetary 
scales. The characteristics and effects described in these 
studies are summarized in the following paragraphs. 

Over northern Asia there is intense, deep cold air 
advection associated with the passage of a mid-tropospheric 
trough. The trough intensifies rapidly over northern Japan. 
The baroclinicity increases with the anticyclogenesis over 
Siberia and cyclogenesis in the strong baroclinic zone off 
the east China coast. As the anticyclone intensifies, 
Subsidence increases which strengthens the east Asian Hadley 
Girculation = winds intensification increases the upper-level 
ageostrophic flow, causing an acceleration of the east Asian 
Subtropical jet near Japan. 

The initiation of the monsoon surge occurs when the upper- 
level trough moves eastward and the pressure gradient along 
the east China coast increases. When the pressure gradient 
force exceeds the Coriolis force, the normal geostrophic 
balance is broken, resulting in Cross-iseobamice fle, scuen 


high to low pressure. This accelerating northeasterly 
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low-level flow toward lower pressure is the monsoon surge. 
The surge reaches the equatorial South China Sea in 12-24 

h. In the equatorial regions, convection intensifies within 
One to two days, sustaining the already enhanced Hadley 
circulation. The upper-level outflow from the convection 
accelerates eastward and westward along the equator as well 
aS northward. The surge ends when the mid-latitude trough/ 
ridge pattern moves far enough east, thus diminishing the 
strong subsidence over the anticyclone. 

The monsoon surge has two stages separated by a few 
hours to one day depending on the location of the observing 
aeons Haneda chang, 1931; Chang et al., 1983). The pres- 
Surge surge is the first stage. It is the leading edge of 
the downward accelerating air with a propagation Speed 
characteristic of internal gravity waves. A frontal passage 
is the second stage which moves with advective speeds and 
is defined by a Sharp decrease in surface dewpoint. 

The cause or nature of monsoon surge forcing has been 
examined by Lim and Chang (1981) and Baker (1983). Linearized 
Shallow-water equations on an equatorial beta-plane were 
used by Lim and Chang (1981) in an analytical study of 
monsoon response to simulated anticyclogenesis. They sug- 
gested that the large features of synoptic scale motions 
in the east ASian winter monsoon could be explained by 
pressure forcing uSing equatorial beta-plane dynamics. No 


consideration was given to orography, a variable mean zonal 


i 


wind, the vertical structure of monsoon surges or the upper 
tropospheric response to the surge. 

Baker (1983) used a global, six-level primitive equation 
model to examine the interaction of a mid-latitude baro- 
clinic wave with topography. A well-developed baroclinic 
wave could initiate a monsoon surge, however, the surges 
were weak with limited extent. Baker (1983) concluded that 
other forcing mechanisms were required to simulate the 
surges observed during the winter monsoon. 

Shaffer et al. (1984) have shown Ehat cold j2urge 2 ie. 
are frequently preceded by rapid growth in the planetary 
waves. It appears that the amplified planetary wave helps 
advect the cold air toward the equator. The purpose of this 
study is to examine the interaction of a planetary wave 
with a growing baroclinic wave. The simulations will be 


carried out with a global, six-level primitive equation 


model. The same baroclinic jet will be used, as was em- 
ployed by Baker (1983). The baroclinic wave, which is unstable 
with the jet, will have wavenumber 8. The planetary wave 


will be represented by wavenumber 4, which is stable with 


the jet. 
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A version of the UCLA potential enstrophy conserving 
general circulation model is used in these experiments. The 
atmosphere is adiabatic, dry and frictionless. The governing 
equations of the model are the primitive equations (listed 
in Appendix A). The model design and dynamic processes 
are described by Arakawa and Lamb (1977,1981). The version 


used for these experiments iS summarized below. 


A. THE HORIZONTAL GRID 
Arakawa's scheme C is used. The prognostic variables-- 
wind (u,v), potential temperature (0) and terrain pressure 


(T)--are staggered horizontally as shown in Figure 2.1. This 


wr8 T.9 





Heg@ceno- Horizontal distribution of variables in 
the UCLA model grid (Arakawa and Lamb, 
Oia) me 


grid arrangement requires more computer time than an unstag- 


gered grid, but it provides a much better structure for the 


is 


shorter waves ariSing from geostrophic adjustment, smaller 
scale heating, nonlinear interaction or topographic effects. 
Grid resolution is 2.5*° longitude by 4.0" latitude.” A ous 
sector of the globe extending from pole to pole is used as 
the computational domain. Cyclic continuity is imposed at 
the east and west boundaries. 

Convergence of the meridians toward the poles normally 
requires a very short time step to ensure computational 
Stability. Instead of shortening the time step, longitu- 
dinal averaging is done on the terms governing gravity wave 
propagation (Arakawa and Lamb, 1977). The averaging is done 
by expanding the zonal pressure gradient and zonal mass flux 
into a Fourier series around each latitude circle. The 
amplitude of each wave component is then reduced by the 
factor S. The gradients of Ehe meteouelocgres amie helo mona. 
smoothed instead of their Fourier expansions. Application 
of the smoothing operator, S, occurs at high latitudes where 
S < l. The smoothing operator is defined as: 

Ax cos 0 


Sn AG Seba ay (2.1) 


(A list of symbol definitions is contained in Appendix B.) 


Bo. “THE VEREPCGAL GRE 
The model equations are written using the vertical sigma 


coordinate: 
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ao Se 
o7 = jr WZ) 


where: 
= Po P,- 


om is the surface pressure, and P, 1s the pressure at the 
top of the model atmosphere (Dy = 200 mb). The variables 


are vertically staggered in six equally spaced sigma layers 























(Figure 2.2). The surface and terrain pressure are defined 
p=200mb o=0O c=0 
a eee Uy Oo = 
o 
a Un = = = 
os 
—_——-— ——-— — uve--— 
o 
— a ae eee ee oe U.v.Q—-— 
o 
—_— —---—-— UVvEe--- 
o 
—_—- ee OU BH 
P=Ps SSS G=0 
PiGuGenz 2. vertical GiStribution of variables in a 
Six-layer sigma coordinate system. Solid 


lines are o-levels; dashed lines represent 
the middle of the layer (Arakawa and 
Eamon oy 7 jc. 


on the lower boundary, o = 1. The vertical velocity, oOo, 


vanishes at the upper and lower boundaries and is diagnostically 


JG 


derived on sigma levels. Wind (u,v) and potential tempera- 


ture (8) are computed at mid-layer. 


C. FINITE DIFFERENCING SCHEMES 

The model uses two time differencing methods--an Euler- 
backward scheme and the centered time differencing (or 
leapfrog) scheme. The model integration proceeds in a 
series of one Eurler-backward time step followed by four 
centered time steps. The Euler-backward scheme eliminates 
the computational mode of the centered scheme and selectively 
damps high frequency waves. The time step for these model 


integrations is 360 s. 
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Mit oben k ENT TIALEZAT ION 


An analytical initialization scheme developed by Hayes 
(1985) is used. In this scheme, the initial mean wind is 
specified analytically using the Bickley jet described by 
Haltiner and Williams (1980). The initial surface pressure 
is derived from the wind using the geostrophic wind eeqacre on, 
and the initial temperature is derived using the thermal 
wind equation. In this research, a stable planetary wave 
and a synoptic scale perturbation are superimposed on a 


horizontally and vertically varying mean current. 


A. INITIALIZATION OF THE MEAN FLOW 
A baroclinic jet is located at 45° N and has a maximum 
velocity of approximately 46 m/s at 200 mb. The initial 


mean surface flow is given by the following equations: 


— — 2 
ug (>) u,, sech (y (-,)) eee 


eee = 0. (SZ) 
It has the same latitudinal variation as the jet with a 
maximum speed of 5 m/s at 45° N decreasing to zero at the 


poles. The v-component of the wind is initially set to 


zero. The mean upper-level wind is given by: 


EY, 


= = &n (ps /p) 
u,, (o,P (9,6) ) ~ , sech* LYCOS °o)) GBBT? (323) 


Vi,(orP(d,6)) = 0 (3.4) 


The wind varies in the vertical as the natural logarithm of 
pressure. 
The mean surface pressure is calculated by integrating 


the geostrophic wind equation: 


- _ p = 

p.(o) = exp[%£n pl(6)) - (a/R) f ful(o)/TS(¢)4¢] (3.5) 
) 0 

where p. (6) = 45° N) = 1013.25 mb and T.(¢) represents the 

mean Surface temperature. Because the latitude structure of 


T (¢) is unknown, p.(9) 1s computed iteratively until the 
adjustment of we to the NACA standard for that Pe is less 
thant Oh 0 Fe mice 

The mean surface temperature is calculated by integrating 


the geostrophic thermal wind equation which gives: 


a _ u_a2o p 2 
T(p,9) = T(P, oq) a R In (p/p,) ; f Sin ¢ sech (d-b,5) de (356) 
0 


The integrals are computed using the Simpson Rule approxima- 


Eon. T (pro) is obtained from the NACA standard atmosphere: 


20 


eee (Py) = 12 Got) 


where: 


The value of T PoC eet S eaken E£rom the NACA standard 


Or 
atmosphere using a Standard surface pressure of 1013.25 mb. 


(Constants are defined in Appendix B.) 


oye INITIALIZATION OF THE HAURWITZ WAVE 

The initialization of the Haurwitz wave is accomplished 
in a Subroutine developed by Dr. T. Rostiond.~ The analyti- 
cal solution for a Rossby wave on the surface of a sphere 
1s described by Haurwitz (1940). 

The Haurwitz wave has a zonal wave number 4 (Figure 
3.1). The maximum meridional velocity is 16 m/s. The u and 


v components are given by the following equations: 


aaa Ow 
0} ADS (S75) 
V' = oe ae (no) 
a cos odd 


The equations of motion are used to derive the perturbation 


pressure: 


SPersonal communication, present affiliation Naval 
Environmental Prediction ReSearch Facility, Monterey, CA. 


2M 


Fagure “3 ie: 





Sea-level pressure field for a Statlonaimy 
Haurwitz wave (3.10). Maximum meridional 
velocity is 16 m/s. Contour interval is 
20 inh: 
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n+l 
n 


n-mt+l m 


P' = C2(w+2)cos(g8t+ma) | ate Syed 





sinop. - (3.10) 
where p, (sin ?) 1S the associated Legendre polynomial and w 
is the mean flow. 

The equations listed above are for a Haurwitz wave 
superimposed on a mean flow. The Bickley jet and the 
Haurwitz wave schemes both have mean flow contributions, and 
one of the mean flows needs to be removed. Instead of using 
the formal analytical approach, the mean u, v and p for 
each J (latitude strip) is subtracted from the Haurwitz solu- 
tion to provide the HaurwitzZ wave contribution. This method 
allows the sea-level pressures to be computed at the same 
time as the winds for the Haurwitz wave and the Bickley jet. 
Otherwise, the winds would have to be computed and then the 
pressures. The error in this approach should not have a 
Significant effect on the results. 

The Haurwitz wave u, v and p contributions are added to 
the Bickley jet and mean surface flow values. The surface 
values are spread to the upper levels using the equations 


described in the previous section. 


om THE PERTURBATION 

The perturbation is a barotropic disturbance consisting 
of a longitudinal sine wave which is used for baroclinic 
instability studies. The maximum latitudinal amplitude in 
the meridional direction occurs at 45° N. The following 


equation describes the perturbation: 


DS 


Es Shia (md) sin” *( 


0 26) ) (3.15 


The hydrostatic equation and the geostrophic wind 


equation are used to derive: 


' —_ a pet — — EOF, . > 2 
p RT RT (A sin (mA)sin (26)) (3. 
(oo. ee eee . 
u = Ne, 36 5 (A sin (mA)2sin(4¢) ) (3 cols 
an a al oo) | 1 ee 
V = af, cos 6 yy = a cos gt BOs (mA) Sin (20) ) (GSalea 
where A = a-l1 m/s. These perturbation quantities are added 


to the mean current and integrated forward in time to 


produce the forecast. 


24 


IV. ANALYSIS OF RESULTS 


The hypothesis is that a cold surge results from the 
interaction of a growing wavenumber 8 cyclone moving east- 
ward and a quaSi-stationary stable forced wavenumber 4. 
Hayes (1985) showed that the basic baroclinic current, which 
is given in Equations (3.1) and (3.3), is unstable for wave 
number 8 but is stable for wavenumber 4. The stationary 
wave could be partially forced by the Himalayas. Three 
experiments were run to verify the hypothesis. Experiment 
I shows the effect of the stable wavenumber 4 on the south- 
ward migration of cold air. The effect of cyclogenesis in 
wave.number 8 is depicted Experiment II. Experiment III 


is a combination of Experiments I and II. 


Ewe EAPERIMENT I 

In this experiment, the model is initialized with the 
stable wavenumber 4 and the baroclinic mean flow described 
mimchiapter IIl. 

The baroclinic jet is confined to the northern hemisphere 
to minimize forcing from the southern hemisphere. The 
model is integrated forward in time to 60 h. Sea-level 
pressure fields for 12, 24, 36, 48 and 60 h are shown in 
Figures 4.1 to 4.5. The potential temperature fields at 
o = 11/12 for the same times are displayed in Figures 4.6 to 


4.10. The corresponding wind vectors are shown in Figures 


2 





Figure 4.1. Sea-level pressure field at 12 h for 
Experiment 1: “COnteumeinter 2 ace omer 
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Sea-~level pressure field at 24 h for 


Figure 4.2. 


Comecure dnterval is 5 mb. 


Experiment I. 
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Figure 4.3. Sea-level pressure field at 36 h for 
| Experiment I. Contour interval is 5 mb. 
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Figure 4.4. Sea-level pressure field at 48 h for 
Experiment I. Contour interval is 5 mb. 
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Figure 4.5. Sea-level pressure field at 60 h for 
Experiment I. Contour interval is 5 mb. 
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Figure 4.6. 





Potential temperature field at 12 h for 
Experiment I. Contour interval is 5° K. 
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Figure 4.7. Potential temperature field at 24 h for 
Experiment I. Contour interval is 5° K. 
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Figure 4.8. Potential temperature field at 36 h for 
PecpeCuimMene et. CONntour interval is 5° K. 
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Figure 4.9. Potential temperature fieldvac 4oeieco. 
Experiment 1... Contour imtervaleises ee 
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Figure 4.10. Potential temperature field at 60 h for 
Peper imene is Contour interval is 5° K. 
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4.11 to 4.15. The 90° sectors are plotted twice toncreduce 
hemispheric displays. 

In the first twelve hours, the effects of the baroclings 
jet on the wave are readily apparent (Figure 4.1). The wave 
has moved westward approximately 5° compared to 10° in the 
southern hemisphere. The westward movement is caused by the 
beta effect. A slight NE-SW tilt in the sea-level pressure 
field is apparent. The absence of a corresponding tilt 
in the southern hemisphere is consistent with the absence of 
the jet. In the northern hemisphere, the tilt is in the 
same sense as the shear on the south side of the jet indi- 
cating that the wave will lose energy to the mean flow. The 
wave should damp. 

The pressure gradient has weakened as the cyclone and 
anticyclone are reduced in amplitude by 5 mb. However, at 
the same time, the amplitude of wave number 4 in the poten- 
tial temperature field (Figure 4.6) has increased. The 
305° K isotherm has moved south to near 10° N. The jet has 
advected that field eastward to a position 45° out of phase 
with the pressure field. Northeasterly winds of approxi- 
mately 13 m/s are observed at 15° N (Figure 4.11). 

By 24 h, the tilt in the sea-level pressure trough has 
increased to produce a stronger pressure gradient in the 
northwest quadrant and a weaker gradient in the northeast 
quadrant of the anticyclone. There has been a corresponding 


Change in the wind field (Figure 4.12). A large area of 


36 


on 
can : d fo NN, A os 
AA mA ee ANE ES 


o 


é 
‘eicce 





60 NK . 
er ies Oo Se ae ¢ mS 
mn. ms eae a so Av AC 7 
Sees Vf 7 ra A ee  / 
SN nae | yi ’ aa A \ aN 7. 
\ ates Se are 
“4 aw, Yes ae ia ek d oP 
{ | a a7G ff dk , | a | \ o~ ¢ / ; 
tee ff Ae eee 
a wf inm od yp \ ae La a 
Das 8- ee De 9° 
Mee ST oN 
-~ eat, ie os ee 
, ~ \ \ / / oS ~~ \ \ alee —~ \ i) i wv” me ; | 
\ : 17 i, SG / 
\ prey , Y Sq ! , ve het) as 1. / 
30° J \--- 7 on ' 7 ! a yh ay aa ise S 
\ / it \ \ rt Athy i \ } 7. 
ae INNA, 7 yar? ine: Poy\ 
~ 


P J 


a 2a mS a SS seca 


90° SP-00"° S 


Figure 4.11. Wind vectors at 12 h for Experiment 
I. 8 mm = 20 m/s. 
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Wind vectors at 24 h for Experi- 
ment I. 7 mm = 20 m/s. 
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Figure 4713. Wind vectors at 36 h for Experi- 
ment I. 7 mm = 20 m/s. 


39 


0 


90 “N 
ee en 
Precick oe SW, ~~. 
1 Sere OL ARN yess Kw. 60° N 


60° Nez 
Cae i“ Shee www y f= S 
/~—0 A ee Few xX 


eee 
20) No : , Pa Ge i ad - 


am) 
, mee j fe ee Y ~. 
\\1 1273 oe mae le a 
/ ‘ a= ‘ ; / — 1. 
ip pf if i a a 2 , f A | ge | 
4.2A N70 CIN TSS A 
a rir ae 4. 735 6 ee 
Gu 4 i \- | j pas ae j i \ : ] ae ' ! Yu j L/ 4 
305 2 a pry, y i \ 7 \ 20 
M44 ean 5 eg! | ‘Ne, | eee / 
NG vsac py asc tsar py cea 


~ ae aay : 


ay be wee FN e 
‘7a is inves ey 
aa ai || || 


ee 


Y ~~. Eee 
90° Se ga0" s 


Figure 4.14. Wind vectors at 48 h for Experi- 
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northeasterly winds is present in the southeast quadrant. 
Pressure centers have continued to weaken. A well-defined 
trough has formed along the equator with wind speeds as 
high as 13 m/s. The tilt and amplitude of the povenenel 
temperature wave have increased. Strong southwesterlies 
ahead of the trough are advecting warm air northward from 
202 IN tO 452 \- 

At 36 h, the anticyclone has intensified by about 5 mb. 
Ridging dominates the pattern particularly south of 30° N. 
The trough has filled as the flow has become more zonal. 

The cyclone has intensified slightly. Weak cyclonic flow is 
evident at the base of the trough (Figure 4.13). Amplifi- 
cation of the potential temperature field has continued with 
strong warm air advection ahead of the trough. Cold air 
advection behind the trough increased south of 20° N. With 
cold air moving south and warm air moving north, V'é' is 
positive, which indicates a baroclinic process. The area 
and speed of northeasterly winds on the eastern side of the 
anticyclone have increased. 

At 48 h, the ridging south of 30° N has expanded east- 
ward while the anticyclone weakened. North of 45° N the 
cyclone has intensified slightly and enlarged. The pressure 
gradient is oriented north-south. Although not evident on 
the pressure field, a weak cyclone has formed at the base of 
the trough (Figures 4.4 and 4.14). The wave structure in the 


potential temperature has amplified further. The 305° K 
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isotherm has moved to within 3° of the equator. North of 
Eee eN,e the £1E€1d is 180° out of phase with that south of 
30° N. Winds are northeasterly on the upwind side of the 
thermal trough. 

In the next twelve hours, rapid development has occurred 
in the cyclone. The central pressure has dropped almost 
10 mb. A secondary anticyclone formed 15° northwest of the 
primary high pressure center. The wave in the potential 
temperature field has continued to amplify. The 305° K 
isotherm has reached 1° S. The field does not reflect the 
surface pressure pattern north of 30° N. Northeasterly 
winds of greater than 15 m/s on the eaStern side of the 
primary anticyclone are observed in the equatorial region 
(Figure 4.15). The secondary anticyclone is not adequately 
reflected in the wind fields indicating that this feature 
is not in geostrophic balance. 

This experiment demonstrates the effects of a baroclinic 
jet on the stable wave number 4. The jet provided enough 
forcing to produce cold temperatures and strong northerly 
winds in the tropics. The 305° K isotherm moved southward 
peom near LO”~ N to south of the equator. At 60 h, the 
secondary pressure centers indicate growth of wave number 8 
which is unstable. This wave is induced by the nonlinear 


interaction of the long wave with itself. 
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By Sela IEMNSIN (SE IE 

In this experiment, the model was initialized with the 
baroclinically unstable mean flow and the wave number 8 
disturbance described in Chapter III. The model was 
again integrated forward in time to 60 h. Sea-level 
pressure fields for 12, 24, 36 and 48 h are shown in Figures 
4.16 through 4.19. The potential temperature fields for 
12, 24, 36 and 48 h are depicted in Figures 4.20 through 4.23. 

The results of this experiment are used to determine 
when to add the stable wave number 4 flow which will advect 
the synoptic scale disturbance southward in Experiment III. 
One requirement for this simulation is to have a sea-level 
pressure difference of at least 25 mb between high and low 
centers to ensure an adequate northwesterly current for 
Surge initiation (Baker, 1983). Another requirement is the 
appearance of a high pressure center representing the 
Siberian anticyclone to the northwest of the low pressure 
Center. 

As shown in Figure 4.16, the disturbance is primarily a 
feature of the northern hemisphere mid- and upper-latitudes 
with a very flat pressure gradient in the tropical region. 
In the first twelve hours a low has formed in the trough 
approximately at 60° N. The trough has a slight NE-SW 
tilt and has moved eastward approximately 10° at 45° N. 

A sea-level pressure difference of greater than 25 mbee tse 


between high and low centers. The potential temperature and 
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Sea-level pressure field at 12 h for 


Figure 4.16. 


GCenmeemws interval 1s 5 mb. 


Experiment II. 
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Sea-level pressure field at 24 h for 


Paigure 4207. 


Conteur interval eis ono 


Experiment II. 
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Sea-level pressure field at 36 h for 


Pigune 4518 . 


Gomteour interval is 5 mb. 


Pxpesiment if. 
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Figure 4.19. Sea-level pressure field at 48 h for 
Experiment II. Contour interval is 5 mb. 
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Figure 4.20. 
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Potential temperature field at 12 h for 
Epemiment ii. “GOntour interval is 5° 
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Figure 4.21. Potential temperature field at 24 h for 
Experiment Ll, Conteum intervals 65. ae 
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Potential temperature field at 36 h for 


Experiment II. 
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Figure 4.23. Potential temperature field at 48 h for 
Experiment il. ~Conmteur Antervaleerse> ace 
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wind fields indicate a baroclinic process with warm air 
being advected northward and cold air being advected 
southward. 

At 24 h, the low centers have moved 5° southeast of 
their previous positions and have intensified. A sharp 
surface-pressure ridge extending north to 80° N separates 
the disturbances. Eastward movement of the trough has 
continued at the ae ee rate. Development of the primary 
disturbance in the wind and potential temperature fields 
parallels that of the pressure field. Some short wave noise 
1S apparent in the wind field. Noise is to be expected 
because no surface pete eional Grssitpatiton, Glffusien Or 
smoothing techniques are included in the model. These 
techniques tend to dissipate surge energy, and the objective 
of this study is the simulation of. intense monsoon surges. 

From 24 h to 36 h, rapid deepening has occurred in the 
cyclone. The central pressure of the low has dropped at 
least 15 mb. The ridge between low centers has weakened, 
but the sea-level pressure difference between high and low 
centers has increased to more than 50 mb. A broad band of 
low pressure appears to dominate the zone between 50° N and 
60° N. Surface warm and cold fronts are clearly evident 
fence temperature field at 36 h (Figure 4.22). The potential 
temperature shows a poleward flux of warm air and a migration 


Siecola air south of 30° N. 


oe, 


The rapid deepening of the surface pressure low con- 
tinued between 36 h and 48 h. The central pressure dropped 
approximately 10 mb. High pressure centers formed near 
40° N, and the sea-level pressure difference between the 
high and low pressure centers iS greater than 70 mb. The 
northern portion of the strong baroclinic zone is advected 
southeastward about 7° latitude. The potential temperature 
field shows some additional migration of cold air south of 
30° N (Figure 4.23) compared to 12 h earlier. At 60h, the 
300° K isotherm follows the 30° N latitude circle, and the 
305° K isotherm 18S a Littles cocurhwer 20 ae 

Analysis of these fields reveals several key features. 
A sea-level pressure difference of greater than 25 mb 
between high and low centers is present from 12 h onward. 
The sea-level pressure difference actually increases 
throughout the run. A een pressure center northwest of the 
low is best depicted at 24 h; it weakens afterwards. A 
strong Siberian high is not observed probably due to the 
presence of a strong jet and lack of topographic and radia- 
tive forcing. The period of most rapid cyclonic develop- 
ment 1S between 24 h and 36 h. The cold air and strong 
northerly winds never penetrate south of 30° N. 

Based on these results, wave number 4 should be added to 
the growing cyclone between 24 h and 36 h. The most rapid 
cyclogenesis is occurring at this point. Adding the larger 
scale wave will further advect the cold air south to help 


initiate the surge - 
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See APERIMENT IIT 

The initialization for this experiment is the same analy- 
tic scheme as Experiment II. However, after 36 h of model 
time, the wave number 4 fields are added. (Note: The 
zOnal mean has been removed from the wave number 4 fields.) 
Sea-level pressure, potential temperature, and wind fields 
are depicted in Figures 4.24 through 4.35. 

The development of the disturbance through 36 h is 
described in Experiment II. The effect of wave number 4 on 
the disturbance is readily apparent when comparing the 36 h 
fields of Experiments II and III. A high pressure center 
of at least 1030 mb forms in the ridge. The primary cyclone 
deepens another 35 mb with the addition of wave number 4. 
This rapid deepening produces a pressure. surge. 

The second stage of the monsoon surge is characterized 
by a sharp drop in surface dewpoint (Baker, 1983). This 
model does not include moisture or its effects, but observa- 
mona studies (Chang et al., 19/79 and Chu, 1978) show that 
this stage can also be associated with a drop in surface 
Eemperature. At 36 h, a trough of cold air extends south- 
ward to approximately 17° N at the base of the high pressure, 
however, the potential temperature gradients are very weak 
south of 30° N. The winds are primarily northerly at 
15 m/s on the eastern side of the anticyclone. 

From 36 h to 48 h the weaker low moved east toward the 


stronger low. The long wave ridge weakened and tilted more 
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Sea-level pressure field at 36 h for 
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Figure 4.26. Sea-level pressure field at 54 h for 
Experiment III. Contour interval is 
PO mb. 
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Figure 4.27. Sea-level pressure field at 66 h for 
Experiment III. Contour interval is 
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Figure 4.28. Potential temperature field at 36 h for 
Experiment III. Contour interval is 
a he 
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Petia 2 9. POtential temperature field at 48 h for 
Experiment III. Contour interval is 
De Ks 
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Figure 4.30. Potential temperature field at 54 h for 
Experiment III. Contour interval is 
So ee 
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Figure 4.32. Wind vectors at 36, fh ference imenr 
III. 4mm = 20 m/s. 
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NE-SW. The small cyclone deepened at least 10 mb. The 
large cyclone deepened approximately 10 mb and moved 
slightly eastward with the mean flow. The tilting and 
weakening of the ridge caused three anticyclonic circulations 
to develop within the ridge particularly toward the north- 
west. At 48 h, winds are northeasterly from 40° N to about 
10° N with a speed of at least 18 m/s. 

The 48 h potential temperature field shows a broad 
tongue of cold air extending southwest near 20° N as if 
spreading over the South China Sea. The furthest extent is 
approximately 15° N. The thermal gradient remains weak south 
of 30° N and especially south of 20° N. A strong baroclinic 
zone 1S depicted behind the cold front. 

Six hours later the secondary low pressure center 1s 
located 5° latitude further southeastward. The primary low 
center remained stationary. The ridge elongated and moved 
westward. This movement resulted in a band of northeasterly 
winds (of approximately 20 m/s) extending southwestward 
from 40° N to at least 10° N (Figure 4.34) as if along the 
eastern coast of China. These northeasterly winds advect 
colder temperatures as far south as 15° N (Figure 4.30). 

The strong baroclinic zone moved approximately 2.5° 
southeastward. 

From 54 h to 66 h the low pressure center moves east- 
ward, and the high pressure center moves westward. The 


pressure gradient south of 30° N decreases. The NE-SW 
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band of northeasterly winds narrowed but continues to 
advect cold temperatures as far south as approximately 

10° N. The strong baroclinic zone continued moving south- 
eastward but at a slower rate. 

Several key observations are made based on this analy- 
Sis. Northeasterly winds move rapidly south to the tropical 
regions in 12h. These winds brought cold mid-latitude 
air to the tropics in 30 h. Compared to Experiment I, the 
northeasterly winds are stronger but do not go quite as 
far south into the equatorial region. The sea-level pressure 
pattern progresses eastward very little because additional 
forcing is necessary to overcome the beta effect on wave 
number 4. The pressure ridge and the cold air tongue are 
better defined than the wind surge. Baroclinic processes 
are present in the results of this experiment south of 


oo N. 
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V. CONCLUSIONS 


The interaction of a stable planetary wave with a 
growing cyclonic disturbance in a mid-latitude baroclinic 
Flow was examined. The purpose of this study was to deter- 
mine the contribution of a quasi-stationary stable wave 
number 4 to the initiation of a monsoon surge as a wave 
number 8 cyclone develops. 

Several circulation features have been identified as 
characteristic of monsoon surges. They include: (1) the 
strengthening of north-northeasterly winds in the South 
China Sea, (2) a sharp drop in surface temperature at obser- 
vation stations in China and Southeast Asia, (3) the shallow 
structure of the surge, (4) the strong baroclinic zone and 
(5) northerly winds in the tropics. In varying degrees, 
each of these characteristics were observed in one or more 
of the experiments. 

In Experiment I, the jet provided enough additional 
forcing to the wave number 4 to produce cold temperatures 
and strong northerly winds in the tropics. The 305° K 
isotherm moved from near 10° N to south of the equator in 
48 h. The winds were initially northerly with a small area 
of northeasterly winds. Throughout the experiment, the winds 
on the east side of the anticyclone became more northeasterly 
and increased in speed in the equatorial region. A strong 


baroclinic zone was observed on the downwind side of the trough. 
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Experiment II showed the evolution of an unstable baro- 
clinic wave which was limited to the mid-latitudes. The 
300° K isotherm moved south approximately 5° latitude in 
48 h, as opposed to 10° and south of the equator in Experi- 
ment I. A strong baroclinic zone, which is prerastorietic 
of a frontal zone, waS evident in the low-level potential 
temperature fields at the mid-latitudes, but it did not 
propagate southward into the tropics. In this experiment 
there were no Significant northerly winds south of 45° N 
and no cold temperatures south of 30° N. 

Several characteristics of monsoon sSurgeS were observed 
in Experiment III in which the wave number 4 and wave number 
8 were combined. North-northeasterly winds reached the 
tropical regions 12 h after the addition of the long wave. 
Cold air moved southward to approximately 10° N within 30 h 
after the addition of the long wave. A strong baroclinic 
zone waS evident in the potential temperature fields. It 
was advected southward to the tropics. Northerly winds were 
observed in the tropics at least as far south as 10° N. 
Experiment III roughly approximated events occurring during 
a cold surge. 

These results suggest that the combination of quasi- 
Stationary stable long wave and a developing mid-latitude 
cyclone initiate a monsoon surge. Additional forcing appears 
to be required to keep wave number 4 stationary. Wave number 


4 also should be modified so that it does not extend to the 


Gan 


equator initially. Moving the jet and the maximum lati- 
tudinal amplitude of the perturbation to 35° N would increase 
the baroclinicity south of 30° N as well as more closely 
model the observed position of the jet during a surge event. 
This study found that a quasi-stationary stable wave number 

4 and an unstable wave number 8 could initiate a monsoon 


surge. 
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APPENDIX A 


MODEL PRIMITIVE EQUATIONS 


Equation of State: 


Hiyiarostatic Equation: 


Horizontal Momentum Equation: 


Ld 





nf) + eknv + vine) - vn 22 = 9 
Soneinuity Equation: 
oi + Yo (IV) a ae = 0 
lst Law of Thermodynamics: 
9(1C ane) ~~ 9 (oC, n 6) 
ce Qn 8) . ——— a 
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APPENDIX B 


List sor sr ibOrS 


Cons came 

radius of the earth 

Cons cane 

specific heat at constant pressure 
Coriolis parameter 

value of Coriolis parameter at 45° N 
gravitational acceleration 

wave number 

pressure (p =p+t+p' + P') 

surface pressure at 45° N, 1013.25 mb 
maximum model pressure height, 125 mb 
Surface pressure 

pressure at tropopause 

gas constant for air 

smoothing operator 

Eemoerature 

Surface temperature at 45° N, 288° K 
Surface temperature 

zonal velocity component (u = u + u’ + U') 
maximum velocity, surface: 5 m/s, upper: 65 M/s 
surface velocity component 


upper level velocity component 
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meridional velocity component (v = v + v' + V') 
terrain height 

specific volume 

frequency 

halfwidth of jet 

dry adiabatic lapse rate, 6.5° K/km 
potential temperature 

longitude 

gal 4159 

normalized pressure 

vertical coordinate 

vertical velocity in sigma coordinate system 
latitude 

latitude of jet 

geopotential, also current terrain height 
angular velocity of the earth 


angular velocity of air relative to the earth 
(mean flow) 


Sirreoam.) funceion 
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